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We study the nuclear enhancement of the transverse momentum imbalance for back-to-back par- 
ticle production in both p+A and e+A collisions. Specifically, we present results for photon+jet and 
photon+hadron production in p+A collisions, di-jet and di-hadron production in e+A collisions, and 
heavy-quark and heavy-meson pair production in both p+A and e+A collisions. We evaluate the 
effect of both initial-state and final-state multiple scattering, which determine the strength of the 
£N) , nuclear-induced transverse momentum imbalance in these processes. We give theoretical predictions 

7—i . for the experimentally relevant kinematic regions in d+Au collisions at RHIC, p+Pb collisions at 

' LHC and e+A collisions at the future EIC and LHeC. 
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^ . I. INTRODUCTION 

Ultra-relativistic nucleus-nucleus (A+A) collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large 
II | Hadron Collider (LHC) have paved the way to studying important properties of a new state-of-matter created in 
Qh such collisions, the quark-gluon plasma (QGP) In A+A reactions both final-state QGP effects and initial-state 
cold nuclear matter effects modify experimental observables relative to the naive binary collision-scaled proton-proton 
i . (p+p) baseline expectation. To disentangle these effects has become a top priority for the heavy ion program. 
Proton-nucleus (p+A) collisions and electron-nucleus (e+A) reactions are, at present, the only tools that provide an 
opportunity to probe experimentally and understand theoretically cold nuclear matter effects without the complication 
^ 1 of final-state interactions in the QGP 0-3 ■ 

Transverse momentum broadening in p+A and e+A collisions is among the most studied cold nuclear matter 
effects. When a fast parton propagates through nuclear matter, it can accumulate additional transverse momentum 
via multiple scattering with the soft partons inside the big nucleus either before or after the hard collision. This 
phenomenon is known as transverse momentum broadening, and it can be probed experimentally through the nuclear 
modification of single inclusive jet or hadron production in e+A collisions 0-0] , the nuclear broadening of Drell-Yan 
di-lepton or W^/Z boson production in p+A collisions [1,0 and the Cronin effect 0, El- One can also study the 
multiple scattering effects through two-particle correlations. For example, the nuclear enhancement of the transverse 
momentum imbalance of di-jet and di-hadron production was shown to be sensitive to both initial-state and final-state 
. £h ' multiple scattering [12[ ■ 

^ Different theoretical approaches have been employed to com put e and describe the nuclear broadening effect. These 

include the dipole approach , the random walk approach [15| , the diagrammatic Glauber multiple scattering [ll| , 

the color glass condensate approach [lj], [l8[ , soft collinear effective theory [l9l-|2]| , and the high- twist power expansion 
approach [5, EiHHHlj]- Some possible connections and relations among different frameworks have been discussed in [2o| . 
Following our previous study [f2j . we will evaluate the nuclear broadening in the formalism that represents multiple 
scattering as contributions to the cross section from higher twist matrix elements in the nuclear state. This framework 
follows a well-established QCD factorization formalism for particle production in p+A collisions (22l. [23 - l26j . and has 
been previously used to describe cold nuclear matter effects, such as energy loss [27rl29| , dynamical shadowing [30l - [32T | 
and broadening effects d, [t| [T3, [H[ • The purpose of our paper is to apply the techniques developed in our previous 
study of di-jet and di-hadron transverse momentum imbalance [12] to new channels, which will be accessible to the 
experiments in the near future. It differs from more generic parton broadening phenomenology in that the color and 
kinematic structures of the hard part are evaluated exactly. In particular, we will study the nuclear enhancement of 
the transverse momentum imbalance for photon+jet (photon+hadron) in p+A, di-jet (di-hadron) in e+A, and heavy- 
quark (heavy-meson) pair production in both p+A and e+A collisions. These two-particle correlation observables can 
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be studied in d+Au collisions at ^/s = 200 GeV at RHIC, the forthcoming p+Pb run at ^/s = 5 TeV at the LHC, 
and at the planned Electron Ion Collider (EIC) and Large Hadron Electron Collider (LHcC). Our paper presents a 
unified formalism to predict theoretically this observable for multiple final-state channels. 

The rest of our manuscript is organized as follows: in Sec. II we study photon+jet (photon+hadron) and heavy- 
quark (heavy-meson) pair production in p+A collisions. The beginning of this section is used to introduce the basic 
definition of the transverse momentum imbalance, which is common to all studied processes. We then take into 
account both initial-state and final-state multiple scattering to compute the nuclear enhancement in the transverse 
momentum imbalance of the produced back-to-back particle pair. In Sec. Ill we study the nuclear broadening of 
the transverse momentum imbalance for di-jet (di-hadron), and heavy-quark (heavy- meson) pair production in c+A 
collisions. In Sec. IV we present our numerical estimate on the nuclear broadening for the relevant kinematics at 
RHIC, LHC and the future EIC and LHeC. We summarize our results in Sec. V. 

II. NUCLEAR ENHANCEMENT OF THE TRANSVERSE MOMENTUM IMBALANCE IN p + A 

COLLISIONS 

The main purpose of this paper is to study the nuclear enhancement of the transverse momentum imbalance for 
back-to-back particle production in both p+A and e+A collisions, h(P') + A(P) -+ h\(p\) + h,2{pi) + X. Here h and 
A are the incoming hadron (or virtual photon) and the nucleus, respectively. h\ and /12 are the produced particles in 
the final state with transverse momenta p±± and P2±- To lowest order in perturbative QCD, the production of these 
two particles arises from hard 2 — > 2 scattering processes. If we denote by z the h and A collision axis, h\ and /12 
are produced approximately back-to- back in the transverse (x, y) plane: pi± rts —p 2 ±. In p+A collisions, however, 
the incoming parton can undergo multiple scattering before the hard collisions. The produced final-state particles 
are also likely to undergo multiple interactions in the big nucleus if they are strongly interacting. Both initial-state 
and final-state multiple scattering lead to acoplanarity, or momentum imbalance of the observed two particles. To 
quantify this effect, we define the transverse momentum imbalance q± as: 

q± =Pi± +P2±, (1) 
and the average transverse momentum squared imbalance 

«y=(/^*plk)/^ (2) 

Here, da/dPS is the differential cross section with dVS representing the relevant phase space, to be defined for each 
process in the corresponding section. For example, in p+A collisions, dVS = dy\dy2dp\ for photon+jet production 
and dVS = dyidy 2 dpi±dp 2 ± for photon+hadron production. 

The enhancement of the transverse momentum imbalance (or nuclear broadening) in h+A (h=p, 7*) collisions 
relative to h+p collisions can be quantified by the difference: 

A(<zi) = (ql) hA - (ql) hp . (3) 

The broadening A(q^) is a result of multiple quark and gluon scattering, and is a direct probe of the nuclear medium 
properties. We now take into account both initial-state and final-state multiple parton interactions to calculate the 
nuclear broadening A(qj_) for photon+jet (photon+hadron) and heavy-quark (heavy- meson) pair production in p+A 
collisions. In the next section we study di-jet (di-hadron) and heavy-quark (heavy-meson) pair production in c+A 
collisions. 



A. Photon+jet (hadron) production in p+A collisions 

1. Photon+jet production 

Consider the following back-to-back photon+jet production process (the photons here and throughout the paper 
are direct photons) in p+A collisions: 

p(P')+A(P)^ 1 (p 1 ) + J(p 2 )+X. (4) 

Here, P', P are the four momentum of the incoming hadron and nucleus (per nucleon) with atomic number A and p\ 
and P2 are the four momentum of the produced final-state photon and jet, respectively. The light-cone components 
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of the final-state particles are given by 
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where 2/1.2 and pi,2± are the rapidities and transverse momenta, respectively. In leading-order collinear factorized 
perturbative QCD, the photon and the jet are produced exactly back-to-back, pu_ = — p2_i_. It is important to realize 
that to this order the jet is idetical to the leading parton. It is only at next-to-leading order (NLO) that the QCD 
structure of the jet starts to play a role in the theoretical description of physics observables. In heavy ion collisions, 
experimental observables that include jets in the final state are presented after subtraction of the uncorrelated soft 
hadronic background. Thus, they would be directly comparable to the results presented in this paper. The differential 
cross section to leading order can be written as |33| 



da 



fa/p( x ').fb/A(x) 



dy 1 dy 2 dp A l 
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(6) 



where ^ a b runs over au possible parton flavors, s = (P' + P) 2 is the center-of-mass energy squared, f a / p and f b /A 
represent the proton and nuclear parton distribution functions, respectively. At this order in perturbation theory, the 
parton momentum fractions x' and x are uniquely related to the rapidities and the jet transverse momentum: 
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Pi_ 



(7) 



H^ b ^ cd (s,t,u) are the partonic cross sections as a function of the usual partonic Mandelstam variables s,i,u. They 




FIG. 1. Leading order Feynman diagrams for photon-jet production, 
are calculated from the Feynman diagrams in Fig. [1] and are given by [33L [34| 
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where N r — 3 is the number of colors. 




FIG. 2. Double scattering diagrams for qq — > 75 (left two) and qg — > yq (right two): (a) and (c) are for initial-state double 
scattering, while (b) and (d) are for final-state double scattering. The blobs represent the tree-level diagrams as shown in 

Fig.m 



In p+A collisions, the energetic incoming parton from the proton can undergo multiple scattering with the soft 
partons inside the nuclear matter before the hard collisions (initial-state multiple scattering). After the hard collisions, 
the leading outgoing parton (opposite the photon) will also undergo multiple interactions in the large nucleus (final- 
state multiple scattering). These interactions lead to an enhancement in the photon+jet transverse momentum 
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imbalance, which can be quantified by A(qj_), as defined in Eq. ([3]). This nuclear broadening A(g^) can be calculated 
in perturbative QCD. A specific method based on double parton scattering has been discussed in detail in Refs. p| 
0, E3, H3| . Our derivation closely follows our previous paper [13] ■ The leading contribution to the nuclear broadening 
comes from the double scattering: either initial-state double scattering, as in Fig. [5] (a) and (c) for the partonic 
channels qq — > and qg — > -fq, respectively; or final-state double scattering, as in Fig. [2] (b) and (d). We calculate 
the contributions from these diagrams in the covariant gauge and obtain the following expression for the nuclear 
broadening of photon+jet production in p+A collisions: 



t.u) 



A (<^> = TT2 7 ■ — r-7^T-7^ „ rr — — . (11) 



^a,o x'x ao— »7<lV ) 1 I 

where T^\{x) — T^) A (x) (or T^? A {x)) are twist-4 quark-gluon (or gluon-gluon) correlation functions associated with 
initial-state multiple scattering, with the following operator definitions 0, 

0111: 

TjiW = J ^-e™P + y- J ^^0( y -- yr )0(-y-)I(p A |F Q +( y2 -)^(O) 7 +^( y -)F+«(yr)bA), (12) 

4W = J ^-e^y I ^^^--tfi-M-Ife-^M^ (13) 

On the other hand, T { J ] A (x) and T ( g Jl(x) are the corresponding twist-4 correlation functions connected with the final- 
state multiple scattering. The y a re given by the same expressions in Eqs. (fT2f and (fl3|) . except for the 0-functions 
that are replaced as follows (9l. Il2j: 

9(y- - yr) ) - V~) %2 )• (14) 

The hard part functions H I ab ^ >ld and d are associated with initial and final-state multiple scattering, respectively; 

and are given by 

( C F H^ ld a=quark 

HL^ d = \ , (15) 

[ C A H l a J b ^ d a^gluon 

{ C FH^ ld d^quark 
H^ jd = < . (16) 

{ C AH^ ld d=gluon 

Here, Cf = (N 2 — 1)/2N C and Ca = -^V c are the quadratic Casimir in the fundamental and adjoint representations of 
SU(3) C , respectively. It is also instructive to recall that the strength of the multiple scattering depends on the color 
representation of the complete scattered parton system [13]. For photon+jet production, since the photon does not 
carry color, the final-state multiple scattering only depends on the color of the jet parton (whether it is a quark or a 
gluon), as can be clearly seen in Eq. (ITB1) . 

2. Photon+hadron production 

For the back-to-back photon+hadron production in p+A collisions, p(P') + A(P) — > j(j>i) + h(p2) + X, the leading 
order differential cross section has the following form: 



da_ _ 2iTa s a em n f ^ fa/ p (x')f b / A (x) vl 



£M*) ' Hab^d(s,t,u), (17) 

XX 

bd 



where the momentum fractions z, x' , and x are given by 

z= V —, x' = ^ (e yi + e V2 ) , x = ^(e-^ + e- y2 ). (18) 
Pi± \/s Vs 

Following our previous paper [13], we can easily generalize the calculation for nuclear broadening A(gj_) in photon+jet 
to photon+hadron production by including the fragmentation function D^/^z), and it is given by 

Mil) = TTT-^ IT _ : „^.,^.., r — — ( 19 ) 
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B. Heavy-quark (heavy-meson) pair production in p+A collision 

1. Heavy-quark pair production 



We now study the heavy-quark pair production, p(P')+A(P) —> Q{pi) + Q{p 2 )+X . At leading order in perturbative 
QCD, the heavy quark Q and anti-quark Q are produced back-to-back through the following partonic channels: 
qq — > QQ, gg — > QQ. Thus, pi± = —p 2 ± and \p\x\ — \P2± \ = P_l- The differential cross section can be written as 



(1(7 



TTa 2 s \ - fa/p(x')fb/A(x) H u 



dyidy 2 dp 
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a , 6 



where the parton momentum fractions x' and x are given by 
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(s,t,u) , 



(20) 
(21) 

p^ + tUq and ttlq the heavy quark mass. The hard part functions H^ b ^Q d , H^ b ^Q d , ^ab->QQ ^ or 

partonic processes a(p a ) + b(p b ) ->• + d(p 2 ), a(pj + &(pb) ->• Q(pi) + d(pa)) o(Po) + KPb) -> Q(Pi) + Qfe) 

relevant to the variable flavor scheme [35| are given in [361 ]. In this paper we work in the fixed flavor scheme (three 
light quarks). Furthermore, if instead of the standard definition of Mandelstam variables one introduces the notation 

s = (Pa+Pb) 2 , i = (p a -pi) 2 - m 2 Q , u = {p b - Pl ) 2 - m 2 Q , 

the two relevant hard part functions can be written compactly as (37j 
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(23) 



The nuclear enhancement of the transverse momentum imbalance A(g 2 ^) in p+A collisions can be easily calculated 
and the final result is given by 
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(24) 



where the hard functions H 1 , and H F , are again associated with initial-state and final-state multiple scat- 
tering and are given by 



fjl _ ttU 
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(25) 
(26) 
(27) 

(28) 



IF IF* 

When the heavy quark mass vtlq — > 0, we recover the published results for H ^ ,^ and H g ' g ^ qq in our previous paper 



2. Heavy-meson pair production 

One can easily extend the above calculation to heavy-meson pair production, such as back-to-back D+D, p(P') 
A(P) — > D(pi) + D(p 2 ) + X. The differential cross section is given by [36] 

da _2-Ka 2 sS -^ f dzi . . ,faJp{x')f b / A {x) v 



dyidy 2 dpi±dp2± 



E / d fD D/Q ( Zl )D m ( Z2 ) H^ QQ (s, i, u), (29) 
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where -D_d/q(zi) and Dq/q(z 2 ) are heavy-meson fragmentation functions, z 2 = z\p 2 ±/pi±, and: 

\/S \/S 



(30) 



with m± = J (pi±/zi) 2 + m,Q. The nuclear broadening A(q^_) is given by 



8^a s \^aMl ^D D/Q ( Zl )D B/Q (z 2 
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(31) 



Including the di-jet and di-hadron production in our previous paper [l2j , we have computed the nuclear broadening 
in the transverse momentum imbalance for all important back-to-back two-particle production channels in p+A 
collisions. We now study two-particle production in e+A collisions in the next section. 



III. NUCLEAR ENHANCEMENT OF THE TRANSVERSE MOMENTUM IMBALANCE IN e + A 

COLLISIONS 

In this section, we use the same approach to calculate the nuclear broadening for di-jet and di-hadron, as well as 
heavy-quark (heavy-meson) pair production in deep inelastic scattering (DIS) of a lepton on a big nucleus, or 7* + A 
collisions. Since the virtual photon does not interact with the soft partons in the nuclear target via the strong force, 
the nuclear broadening only comes from the final-state multiple scattering. 



A. Di-jet (di-hadron) production in DIS 

1. Di-jet production 



Nuclear enhancement of the di-jet transverse momentum imbalance in photo-production has been calculated in 



Ref. [22j. Here we will generalize result to di-jet production in DIS, 

7 *(P r ) + A(P) -> ,h{ Pl ) + J 2 (p 2 ) + X 



(32) 



where the incoming virtual-photon 7* carries momentum P 7 » with virtuality P 2 * = — Q 2 . P is the four momentum of 
the target nucleus, p\ and p 2 are the momenta of final-state jets J\ and J 2 , respectively. We will work in the center 
of mass frame of 7*+ A, in which the light-cone components of the incoming particles are 
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(33) 



At leading order, the two jets have opposite transverse momentum but with the same magnitude p± , and the differential 
cross section can be written as 
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Vs/p± 

with xb = Q 2 /2P ■ P 7 » = Q 2 /(s + Q 2 ). y\ is the rapidity of the first jet J 1; and the rapidity y 2 of the second jet J 2 
is related to y\ as follows: 
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The hard function H^ tb ^ cd is given by 
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where s, t, and u is denned as 

s = (P^+xP) 2 , i={P r - Pl ) 2 , u={xP- Pl ) 2 . (39) 

Since the virtual photon does not interact strongly with the nucleus, the broadening A(<7j_) is only sensitive to 
final-state multiple scattering and is given by 

W-V E4/«.4WHy. l ^ a (s,<,i,0 2 ) 

The hard function H^, b ^ cd can be written as 

C C F H^ b ^ cd b=quark 

<^cd = { ■ (41) 

I C A H^ cd b=gluon 

These expressions suggest that even though H^, b ^ cd are the hard functions associated with final-state multiple 
scattering, the strength of the broadening depends on the color representation of the initial-state parton b: the color 
factor Cf (or Ca ) corresponds to the incoming quark (or gluon). This is not surprising since the rescattering effects are 
only sensitive to the total color of the final two-parton composite state, which is equal to the color of the initial parton 
b (as 7* carries no color). It is easy to show that by setting Q 2 — > 0, we recover the result for di-jet photo-production 
derived in Rcf. 221. 



2. Di-hadron production 

For back- to-back hadron pair production, 7*(f 7 *) + A(P) — > hi(p\) + /i2(f>2) +X, the differential cross section can 
be written as 

Z^Dhi/cWDha/dfo)— H^ b ^ cd {s,t,u,Q z ), (42) 



,2 



dyidy 2 dp 1 j_dp 2 j_ ( s + Q 2 ) , 
where the momentum fractions Zi, z<i, and x are given by 



;|= ^ (e ^ +e ^), Z2 = ?^( e in +e * ) , x = x B + I e ll+l? . (43) 
's y/s s + Q z e yi + e y2 



The nuclear broadening A(q 2 _) has the following form: 



B. Heavy-quark (heavy-meson) pair production in DIS 

1. Heavy-quark pair production 

We now study the heavy-quark pair production, 7*(P 7 *) + A(P) — ► Q(pi) + Q(p2) + X. At leading order, there is 
only one partonic channel that contributes, 7*+.<?— >Q + Q [38j |. The differential cross section can be written as 



da _ ira s a em 1 f g /A( x ) „v , , , ,i 

dy x dp\ ~ ( s + Q2f I - n^ e m x Ta^QQ 



H^ a ^ o0 (s,i,u,Q 2 ,m 2 ), (45) 



where y\ is the rapidity of the heavy quark Q, while the heavy anti-quark Q has rapidity yi — ]n(^/s/mx — e Vl ). The 
momentum fraction x is given by 
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with the transverse mass 
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ttiq. The hard part functions H » 9 _^qq has the following form, 
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Once again s, t, and u differ slightly form the standard Mandelstam variables and are given by 

s = (Rj. +xP) 2 , t= (P 7 . -pi) 2 -m 2 Q , u={xP-p 1 f 
The nuclear enhancement of the transverse momentum imbalance can be written as 
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where H F * g ^.QQ = t qq w ith Ca reflective of the color representation of the QQ system. 



2. Heavy-meson pair production 



(47) 



(48) 



(49) 



One can also easily generalize the above calculation to heavy-meson pair production, such as back-to-back D + D 
production, 7*(P 7 „) + A(P) — > D(pi) + D(p 2 ) + X. The differential cross section is given by 
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2na s a e 



dy 1 dy2dpi±dp2± ( s + Q 2 ) [ 



■D d /q(zi)D D/ q(z 2 



{ fg/A{x) H u 



i*g->QQ 



{s,i,u, Q 2 ) , 



where the momentum fractions z\ and z 2 in heavy-meson fragmentation functions are given by 

Pi± 



(e yi +e V2 ) 



z 2 = ^(e y i+ e y*) 

rJs 



Here the parameter r depends on the heavy quark mass tuq: 



The nuclear broadening for this final state is given by 



Mil) = 
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(50) 



(51) 



(52) 



(53) 



We have now completed all evaluations of the nuclear enhancement of the transverse momentum imbalance in back- 
to-back two-particle production in both p+A and e+A collisions. We will use these expressions to present predictions 
for A(q^) relevant to future experimental measurements in the next section. 



IV. NUMERICAL RESULTS 



In this section we present phenomenological applications of our results. Specifically, we give theoretical predic- 
tions for the nuclear broadening (or enhancement of the transverse momentum imbalance) for back-to-back particle 
production in d+Au collisions at RHIC, for the forthcoming p+Pb collisions at LHC, and for the e+A collisions at 
the planned EIC and LHeC. The only new unknown ingredient in our calculation is the twist-4 quark-gluon and 
gluon-gluon correlation functions. Following [12l |30L l3l|. we parametrize them as follows: 

^ = T ™> ] = e ( Al/3 - *) f ^ A(x) > (54) 

where .f q , g /A{ x ) is the standard leading- twist parton distribution function for quarks or gluons, respectively. In Eq. (|54[) 
£ 2 = 0.12 GeV 2 represents a characteristic scale of parton multiple scattering and was extracted from deep inelastic 
scattering data [3fJ. The definition in Eq. ([54]) is such that £ 2 (A 1 / 3 — l) can be thought of a dynamical quark mass 
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generated in the background soft gluon field of the nucleus in minimum-bias reactions [39(. Our implementation has 
also been successful in describing the nuclear suppression of single inclusive hadron production [SJ H(| and the di- 
hadron transverse momentum imbalance and correlations in d+Au collisions at forward rapidities at RHIC y/s = 200 
GeV [Hj]. 

As demonstrated in Ref. [301 ] . by approximately decomposing a nuclear state into a product of nucleon states the 
parameter £ 2 can be expressed in terms of an averaged gluon field strength squared (F +a F + a ) ~ lim^^o x fg/A{ x )i 
the soft-gluon number density. In this picture, £ 2 represents the strength of the multiple scattering, thus proportional 
to the number of the soft gluons in the nuclear medium. While at tree level f 2 is a fixed number, higher order 
corrections may provide energy dependence to this parameter (for example, at LHC energies £ 2 can be larger in 
comparison to the one used RHIC energies) . We study this possibility phenomenologically by presenting results for a 
range of £ 2 from 0.12 GeV 2 to 0.20 GeV 2 at the LHC (yellow band). For the upper limit we take guidance from the 
growth of the inelastic scattering cross section, from erf™ = 42 mb at y/s = 200 GeV to of^ — 70 mb at y/s = 5 
TeV [i^, The energy dependence of £ 2 (if any) can be tested and further constrained once experimental data at 
different center-of-mass energies become available lj. To estimate the transverse momentum broadening for the new 
channels derived in the last two sections, we use the CTEQ6L for nucleon parton distribution functions (44j |. and 
EPS08 parametrization for nuclear parton distribution functions (nPDFs) [45|. Since the nuclear broadening A(q^_) 
is a ratio of cross sections as in Eq. ((3]), it has very weak dependence on the parametrization for nPDFs. 




FIG. 3. Nuclear broadening A(q5_) for back-to-back heavy-quark pairs (left) and di-jets (right) in p+A collisions as a function of 
binary collision number JV C oii. We choose rapidities y 1 = y 2 = 2 for the y/s = 5 TeV LHC p+Pb collisions and y\ = j/2 = 1 for the 
y/s = 200 GeV RHIC d+Au collisions. For the LHC, the jet transverse momentum is integrated over 30 GeV < p± < 40 GeV, 
while for RHIC the jet transverse momentum is integrated over 15 GeV < p± < 25 GeV. The yellow band is for LHC kinematics, 
with the band representing a variation of £ 2 parameter from 0.12 GeV 2 to 0.20 GeV 2 . The red solid curve is for RHIC kinematics 
with £ 2 = 0.12 GeV 2 . 

In Fig. [3] (left) we plot the nuclear enhancement of the the transverse momentum imbalance A(g^_) for a back-to- 
back heavy-quark pair in p+A collisions as a function of the binary collision number N co \\. To take into account the 
centrality dependence, we have replaced (A 1 / 3 — 1) by (A 1 / 3 — l)(A co n(6))/(A col i(6 min . bias )). We choose the rapidities 
of both the heavy quark and the anti-quark to be y\ = yi — 2, and integrate over the transverse momentum 30 GeV 
< p ± < 40 GeV for LHC p+Pb run at y/s = 5 TeV, and y x = y 2 = 1 and 15 GeV < p± < 25 GeV for RHIC 
d+Au run at y/s = 200 GeV. The yellow band is for LHC kinematics, with the band representing a variation of 
the £ 2 parameter from 0.12 GeV 2 to 0.20 GeV 2 as we discussed above. The red solid curve is for RHIC kinematics 
with £ 2 = 0.12 GeV 2 . For comparison, we also plot in Fig. [3] (right) the nuclear broadening for di-jet production 
using the same kinematics. We find that the nuclear broadening A(g 2 _) is slightly stronger for di-jet production when 
compared to heavy-quark pair production. This is because there are more partonic channels that contribute to the 
di-jet. In particular, the gg — > gg channel is the most important and generates the largest broadening [12}. However, 
this channel does not contribute to heavy-quark pair production. Since the gg —¥ gg channel becomes more important 
at larger center-of-mass energies, we expect the difference in A(q 2 ± ) between heavy-quark pair production and di-jet 



1 Recent comparison of theoretical predictions for the cross section modification in p+Pb reactions [4^1 to new ALICE preliminary data 
|4:jf does not favor a growth of £ 2 . 
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FIG. 4. Left panel: nuclear broadening A{q±) f° r back-to-back photon+jet production in p+A collisions as a function of the 
binary collision number iVcoli- We choose rapidities y\ — y% = 2 for the y/s = 5 TeV LHC p+Pb collisions and y\ — j/2 = 1 
for the y/s = 200 GeV RHIC d+Au collisions. For the LHC, the jet transverse momentum is integrated over 30 GeV < p± < 
40 GeV, while for RHIC the jet transverse momentum is integrated over 15 GeV < p± < 25 GeV. Right panel: same as left 
plot, but now for back-to-back photon+hadron production. For LHC kinematics, we integrate over 10 < p 1 ±_ < 20 GeV and 
5 < Ph± < 10 GeV. For RHIC kinematics, we integrate over 5 < p y ± < 15 GeV and 5 < pu± < 10 GeV. The yellow band is for 
LHC kinematics, with the band representing a variation of f; 2 parameter from 0.12 GeV 2 to 0.20 GeV 2 . The red solid curve is 
for RHIC kinematics with £ 2 = 0.12 GeV 2 . 



production to become slightly larger in going form RHIC to the LHC. This can be seen clearly in our predictions in 
Fig.H 

In Fig. S] (left) we plot A{qj_) for the photon+jet final state in p+A collision as a function of N co \\. Compared to 
the heavy quark or di-jet production in Fig. [31 the nuclear broadening is much smaller. This is because there is much 
stronger final-state multiple scattering in di-jet or heavy quark pair production: both outgoing partons interact with 
the cold nuclear matter. The difference in A(<jQ.) between photon+jet and di-jet (or heavy-quark pair) production is a 
direct prediction of our approach. A comparative experimental study facilitated by future experimental measurements 
will be a very useful test of our formalism. In Fig. [4] (right), we also provide a plot for back-to-back photon+hadron 
production in p+A collisions as a function of N co \\. We find that the magnitude of the nuclear enhancement in 
the transverse momentum imbalance is very similar between photon+jet and photon+hadron production. This is 
expected since the nuclear broadening as defined in Eqs. ([2]) and ([3]) is a ratio of cross sections and should not be 
affected much by the fragmentaiton function (the hadronization process). 

Nuclear-enhanced transverse momentum imbalance in di-jet and heavy-quark pair production in DIS can provide 
a clean measurement of final-state rescattering effect. In Fig. [5] we plot A(q 2 _) for both di-jet and heavy-quark pair 
production in 7* + A collisions as a function of the atomic number A. In the left panel we present predictions for 
the relevant kinematics region for the planned future EIC. We choose the 7* + A center-of-mass energy ^=100 GeV 
and the photon virtuality Q 2 = 16 GeV 2 . Note that the 7*- A center-of-mass energy corresponds to the standard DIS 
kinematic variable W: s — (P 7 » + P) 2 = W 2 . We choose the jet rapidity to be y\ — 2 and integrate the jet transverse 
momentum over 5 GeV < p± < 10 GeV. In the right panel we present predictions for the relevant kinematics region 
for the planned future LHeC. We choose the 7* + A center-of-mass energy -^=800 GeV and the photon virtuality 
Q 2 = 25 GeV 2 , and integrate the jet transverse momentum over 10 GeV < p± < 20 GeV. Our numerical results show 
that contrary to the case in p+A collisions, the growth of the transverse momentum imbalance for heavy-quark pair 
production is larger than that for di-jet production because of the different color factors of the rescattering final-state. 
In the heavy-quark case, only one process is present at leading order and the initial-state parton is a gluon. This 
gives rise to a color factor Ca- On the other hand, there are two processes that contribute to the di-jet transverse 
momentum imbalance and they give rise to a combination of color factors: C'f and Ca- Consequently, the nuclear 
broadening is smaller than the one observed in heavy-quark pair production. 
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FIG. 5. Transverse momentum imbalance increase A(g5_} in 7* + A collision as a function of the atomic number A at 7*-A 
center-of-mass energy. Left panel: y/s = 100 GeV and photon virtuality Q 2 =16 GeV 2 for typical EIC kinematics. We choose the 
jet rapidity yi = 2 and have integrated the jet transverse momentum over 5 GeV < p± < 10 GeV. Right panel: \fs = 800 GeV 
and photon virtuality Q 2 =25 GeV 2 for typical LHeC kinematics [3]. We choose the jet rapidity y\ = 2 and have integrated the 
jet transverse momentum over 10 GeV < p± < 20 GeV. Note that the 7*-A center of mass energy corresponds to the standard 
DIS kinematic variable W: s = (P 7 . + P) 2 = W 2 . The solid curve is for heavy-quark pair production, while the dashed curve 
is for di-jet production. 



V. SUMMARY 



Within a high-twist approach to parton interactions in cold nuclear matter we studied the nuclear enhancement 
of the transverse momentum imbalance for photon+jet and photon+hadron production in p+A collisions, di-jet and 
di-hadron production in e+A collisions, and heavy-quark (heavy-meson) pair production in both p+A and c+A 
collisions. By taking into account both initial-state and final-state multiple scattering, we derived results to lowest 
order in perturbative QCD for the increase in the transverse momentum imbalance of two particle-production for these 
channels. We presented numerical predictions for the kinematic regions relevant to d+Au collision at RHIC, p+Pb 
collisions at LHC, and e+A collisions at a future EIC and LHeC. We found that the nuclear broadening in photon+jet 
(photon+hadron) production is much smaller than the one in di-jet (or heavy-quark pair) production, due to weaker 
final-state multiple scattering. It is also interesting to notice that in p+A collisions the di-jet accumulates more 
nuclear-induced transverse momentum imbalance than the heavy-quark pair production, while in e+A collisions it is 
the other way around. The difference in the nuclear broadening among the various channels is a direct predictions of 
our approach. We emphasize that a comparative study of the transverse momentum imbalance of back-to-back particle 
production, facilitated by future experimental measurements, will be a valuable probe for the multiple scattering effect 
in cold nuclear matter and a test of our theoretical formalism. 
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